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(54) Ultrasound imaging using coded excitation on transmit and selective filtering on receive 



(57) A method and an apparatus for selectively per- 
forming contrast harmonic imaging, tissue harmonic im- 
aging and B-mode flow imaging with an ultrasound im- 
aging system utilizes phase-coded excitation on trans- 
mit (24', 36,38) and selective firing-to-firing, i.e., "slow- 
time", filtering (46) on receive. "Slow-time" filtering in 



combination with transmit phases which change over 
the set of transmit firings results in different effective 
"slow-time" filters corresponding to the different har- 
monic modes within the reflected signal. The transmit 
phases and the "slow-time" filter weightings are de- 
signed to selectively enhance the desired modes while 
substantially suppressing other modes. 
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Description 

[0001] This invention generally relates to ultrasound 
imaging of the human anatomy for the purpose of med- 
ical diagnosis. In particular, the invention relates to 
methods and apparatus for imaging tissue harmonics 
and methods and apparatus for imaging fluid flow with 
or without contrast agents. 

[0002] Conventional ultrasound scanners create two- 
dimensional B-mode images of tissue in which the 
brightness of a pixel is based on the intensity of the echo 
return. Conventional B-mode images are formed from a 
combination of fundamental and harmonic signal com- 
ponents, the former being direct echoes of the transmit- 
ted pulse and the latter being generated in a nonlinear 
medium such as tissue from finite-amplitude ultrasound 
propagation. In certain instances, e.g., obese patients, 
ultrasound images can be improved by suppressing the 
fundamental and emphasizing the harmonic signal com- 
ponents. 

[0003] Tissue harmonic imaging was proposed in an 
article by Averkiou et a!, entitled, "A new imaging tech- 
nique based on the nonlinear properties of tissues," 
Proc. 1997 IEEE Ultrasonic Symp. Propagation of 
sound beams in biological tissues is known to be non- 
linear, giving rise to the generation of harmonics. In tis- 
sue harmonic imaging, energy is transmitted at a funda- 
mental frequency f 0 and an image is formed with energy 
at the second harmonic 2/ 0 . Some of the characteristics 
of the nonlinearly generated second harmonic beams 
are a narrower beam, lower sidelobes than the funda- 
mental and beam formation in a cumulative process, i. 
e . the second harmonic continually draws energy from 
the fundamental during propagation. These character- 
istics contribute to axial resolution improvements, re- 
duction of multiple reflections due to tough windows, 
and clutter reduction due to inhomogeneities in the tis- 
sue and skin layers. 

[0004] Conventional ultrasound imaging systems al- 
so have a so-called "color flow" mode in which the flow 
of blood or movement of tissue can be imaged. Conven- 
tional ultrasound flow imaging methods use either the 
Doppler principle or a time-domain cross-correlation 
method to estimate the average flow velocity, which is 
then displayed in color overlaid on a B-mode image. 
[0005] Measurement of blood flow in the heart and 
vessels using the Doppler effect is well known. The fre- 
quency shilt ol backscaltered ultrasound waves may be 
used to measure the velocity of the backscatterers from 
tissue or blood. The change or shift in backscattered fre- 
quency increases when blood flows toward the trans- 
ducer and decreases when blood flows away from the 
transducer. The Doppler shift may be processed to es- 
timate the average flow velocity, which is displayed us- 
ing different colors to represent speed and direction of 
flow. The color flow velocity mode displays hundreds of 
adjacent sample volumes simultaneously, all color-cod- 
ed to represent each sample volume's velocity. 



[0006] In accordance with a known imaging system, 
the color flow mode employs multiple transmit firings for 
each focal point. Operating on a packet of as many as 
16 transmits, a high-pass wall filter rejects echoes from 

s slow-moving tissue or vessel walls to reduce the signal 
dynamic range for subsequent flow processing, using 
the Kasai autocorrelation algorithm or a cross-correla- 
tion algorithm to estimate the average flow velocity. 
[0007] Although quantitative velocity information may 

io be obtained in conventional color-flow imaging, the abil- 
ity to see physical flow is limited by its clutter rejection 
capability, resolution, frame rate, and axial-only flow 
sensitivity. 

[0008] Digital subtraction methods have been previ- 
15 ously proposed to image moving reflectors in B-mode 
imaging (see Ishihara et ai, "Path Lines in Blood Flow 
Using High-Speed Digital Subtraction Echography," 
Proc. 1992 IEEE Ultrason. Symp., pp. 1 277-1 280 : and 
Ishihara et ai, "High-Speed Digital Subtraction Echog- 
20 raphy: Principle and Preliminary Application to Arterio- 
sclerosis, Arrhythmia and Blood Flow visualization," 
Proc. 1990 IEEE Ultrason. Symp., pp. 1473-1476). 
However, these methods use frame-to-frame subtrac- 
tion, which results in a wall filter having an extremely low 
25 cutoff frequency. The low cutoff frequency is due to the 
long time delay between adjacent frames, which does 
not adequately suppress signals from slow-moving tis- 
sue or vessel walls. 

[0009] U.S. Patent No. 5,632,277 to Chapman et at. 

30 discloses a nonlinear imaging system using phase in- 
version subtraction. The Chapman patent uses "first and 
second ultrasound pulses that are alternatively transmit- 
ted into the specimen being imaged," and mentions the 
particular embodiment of transmitting and summing on 

35 receive two pulses that differ by 180°. 

[0010] Contrast agents have been developed for 
medical ultrasound to aid in diagnosis of traditionally dif- 
ficult -to- image vascular anatomy. For example, the use 
of contrast agents is discussed by de Jong etal. in "Prin- 

40 ciples and Recent Developments in Ultrasound Con- 
trast Agents," Ultrasonics, Vol. 29, pp. 324-380 (1991). 
The agents, which are typically microbubbles whose di- 
ameter is in the range of 1-10 micrometers, are injected 
into the blood stream. Since the backscatter signal of 

45 the microbubbles is much larger than that of blood cells, 
the microbubbles are used as markers to allow imaging 
of blood flow. One method to further isolate echoes from 
these agents is to use the (sub)harmonic components 
of the contrast echo, which is much larger than the har- 

50 monic components of the surrounding tissue without 
contrast agent. [See, e.g., Newhouse etal, "Second 
Harmonic Doppler Ultrasound Blood Perfusion Meas- 
urement," Proc. 1992 IEEE Ultrason. Symp., pp. 
1175-1177; and Burns, et at., "Harmonic Power Mode 

55 Doppler Using Microbubble Contrast Agents: An Im- 
proved Method for Small Vessel Flow Imaging," Proc. 
1994 IEEE Ultrason. Symp., pp. 1547-1550 ] 
[0011] U.S. Patent No. 5,706,81 9 to Hwang et ai. dis- 
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closes a method and an apparatus for ultrasonic imag- 
ing using harmonic contrast agents, e.g., gas-filled mi- 
crobubbles. Ultrasonic pulses of opposite polarity are 
transmitted in successive firings. The respective echo 
signals are summed to extract the harmonic response 
attributable to the injected contrast agents. 
[0012] There is a need for a method of contrast har- 
monic imaging in which fundamental or second harmon- 
ic signals from contrast flow can be visualized with sup- 
pressed background tissue signals and little motion 
flash artifacts. This requires the imaging system to have 
high dynamic range, the ability to reject clutter from sta- 
tionary or slow moving tissue and vessel walls, high res- 
olution, high frame rate, and flow sensitivity in ail direc- 
tions. A need also exists for a method of tissue harmonic 
imaging in which harmonic signals generated by nonlin- 
ear propagation in tissue are visualized. In addition, a 
need exists for a method ol visualizing fundamental sig- 
nals from blood flow (without contrast agents) in B mode 
with minimal motion flash artifacts. A further need exists 
for a programmable ultrasound imaging system capable 
of selectively performing contrast harmonic imaging, tis- 
sue harmonic imaging and B-mode flow imaging. 
[0013] A method and apparatus are provided for se- 
lectively performing contrast harmonic imaging, tissue 
harmonic imaging and B-mode flow imaging wherein the 
preferred embodiments use phase-coded excitation on 
transmit and selective firing-to-firing, i.e., "slow-time 0 , 
filtering on receive. "Slow-time" filtering in combination 
with transmit phases which change over the set of trans- 
mit firings results in different effective "slow-time" filters 
corresponding to the different modes within the reflected 
signal. The transmit phases and the "slow-time" filter 
weightings are designed to selectively enhance the de- 
sired modes while suppressing others. In particular, a 
sequence of broadband pulses with different phases 
(and possibly different amplitudes) are transmitted to a 
transmit focal position over multiple firings, and the set 
of received beamformed signals are multiplied with a set 
of (possibly complex) scalar weightings before summing 
together that set of weighted beamformed signals for 
subsequent processing to form one image scan line. A 
complete image is formed by repeating this procedure 
for multiple transmit focal positions across the region of 
interest. 

[0014] In accordance with a preferred embodiment of 
the invention, the "slow-time" filter is embodied as a fi- 
nite impulse response (FIR) filter which receives a first 
set of filter coefficients for filtering the receive signal pro- 
duced as a result of a first phase-encoded transmit fir- 
ing, and which receives a second set of filter coefficients 
for filtering the receive signal produced as a result of a 
second phase-encoded transmit firing. The first set of 
filter coefficients is formed by multiplying each of a pre- 
determined set of filter coefficients by a first scalar 
weighting; the second set of filter coefficients is formed 
by multiplying each of the predetermined set of filter co- 
efficients by a second scalar weighting. The transmit 



phases and the "slow-time" scalar weightings are pro- 
grammable as a function of the three different applica- 
tions, namely, contrast harmonic imaging, tissue har- 
monic imaging and B-mode flow imaging. 

5 [0015] In contrast harmonic imaging, the sonographer 
wishes to see fundamental or second harmonic signals 
from contrast flow with suppressed background tissue 
signals and little motion flash artifacts. This can be 
achieved in the following ways: (1 ) by high-pass filtering 

to the second harmonic and suppressing a substantial 
fraction of the fundamental signal, which results in good 
suppression of background tissue signals; (2) by high- 
pass filtering the fundamental and second harmonic sig- 
nals, which results in good background suppression but 

is larger motion flash artifacts; or (3) by high-pass filtering 
or suppressing the fundamental and all-pass filtering the 
second harmonic, which results in more tissue back- 
ground (from the second harmonic) but shows harmonic 
signals from even the slowest -moving contrast agents. 

20 [0016] In tissue harmonic imaging, the goal is to see 
harmonic signals (in particular, the second harmonic) 
generated by nonlinear propagation in tissue. This is 
achieved by suppressing a substantial fraction of the 
fundamental signal and passing a substantial fraction of 

25 the second harmonic signal. 

[0017] Finally, the goal in B-mode flow imaging is to 
visualize fundamental signals from blood flow (without 
contrast agents) with minimal motion flash artifacts. This 
is accomplished by high-pass filtering the fundamental 

30 and all-pass filtering the second harmonic, which 
smoothes out flash artifacts. 

[0018] The "slow-time" filtering is preferably per- 
formed by an FIR filter with B-mode image feed-through. 
The "slow-time" or firing-to-firing filtering permits a long- 

35 er FIR filter for better clutter suppression, while increas- 
ing the cutoff frequency to a useful range. 
[0019] Embodiments of the invention will now be de- 
scribed, by way of example, with reference to the ac- 
companying drawings, in which: 

40 [0020] FIG. 1 is a block diagram of a conventional ul- 
trasound imaging system. 

[0021] FIG. 2 is a block diagram of portions of an ul- 
trasound imaging system in accordance with a preferred 
embodiment of the invention 

45 [0022] FIG. 3 is a flowchart showing "slow-time" filter- 
ing with B-mode feed-through in accordance with anoth- 
er preferred embodiment of the invention. 
[0023] FIGS. 4-11 are graphs showing the filler re- 
sponse (as a function of slow-time normalized frequen- 

50 cy) to the fundamental mode (solid lines), the second 
harmonic (dashed lines) and the second subharmonic 
(dotted lines). The transmit phases and the "slow-time" 
filter weightings are as follows: FIG. 4) [0°, 180°, 0°, 
180°] and [0.4, 1, 1,0.4]; FIG. 5) [0°, 90°, 0°, 180°] and 

55 [0.4, 1,1, 0.4] (with "slow-time" filter phases [0°, 90°, 0°, 
0°]); FIG. 6) [0°, 180° : 180°, 0°] and [0.4, 1,-1, -0.4]; 
FIG. 7)[0°, 180°] and [1, 1]; FIG 8) [180°, 0°, 180°] and 
[0.5, 1,0.5]; FIG. 9) [0°, 0°, 180°, 180°] and [1, 1,1.1]; 
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FIG. 10) [0°, 180°, 180°, 0°] and [1, 1, 1, 1]; and FIG. 
11) [0°, 180°, 0°, 180°] and [1, 1, -1, -1]. 
[0024] The present invention can be incorporated in 
an ultrasonic imaging system of th type depicted in 
FIG. 1 . This imaging system comprises a transducer ar- 
ray 10 comprising a plurality of separately driven trans- 
ducer elements 12, each of which produces a burst of 
ultrasonic energy when energized by a pulsed wave- 
form produced by a transmitter 14. The ultrasonic ener- 
gy reflected back to transducer array 10 from the object 
under study is converted to an electrical signal by each 
receiving transducer element 1 2 and applied separately 
to a receiver 16 through a set of transmit/receive (T/R) 
switches 18. T/R switches 1 8 are typically diodes which 
protect the receive electronics from the high voltages 
generated by the transmit electronics. The transmit sig- 
nal causes the diodes to shut off or limit the signal to the 
receiver. Transmitter 1 4 and receiver 1 6 are operated 
under control of a master controller 20 responsive to 
commands by a human operator. A complete scan is 
performed by acquiring a series of echoes in which 
transmitter 14 is gated ON momentarily to energize 
each transducer element 12, and the subsequent echo 
signals produced by each transducer element 1 2 are ap- 
plied to receiver 16. A channel may begin reception 
while another channel is still transmitting. Receiver 16 
combines the separate echo signals from each trans- 
ducer element to produce a single echo signal which is 
used to produce a line in an image on a display monitor 
22. 

[0025] Under the direction of master controller 20, 
transmitter 14 drives transducer array 10 such that the 
ultrasonic energy is transmitted as a directed focused 
beam. To accomplish this, respective time delays are 
imparted to a plurality of pulsers 24 by a transmit beam- 
former 26. Master controller 20 determines the condi- 
tions under which the acoustic pulses will be transmit- 
ted. With this information, transmit beamformer 26 de- 
termines the timing and amplitudes of each of the trans- 
mit pulses to be generated by pulsers 24. The ampli- 
tudes of each transmit pulse are generated by an apo- 
dization generation circuit 36, which may be a high-volt- 
age controller that sets the power supply voltage to each 
pulser. Pulsers 24 in turn send the transmit pulses to 
each of elements 12 of transducer array 10 via T/R 
switches 1 8, which protect time-gain control (TGC) am- 
plifiers 28 from the high voltages which may exist at the 
transducer array. By appropriately adjusting the transmit 
focus time delays in a conventional manner and also ad- 
justing the apodization weightings, an ultrasonic beam 
can be directed and focused to form a transmit beam. 
[0026] The echo signals produced by each burst of 
ultrasonic energy reflect from objects located at succes- 
sive ranges along each transmit beam. The echo signals 
are sensed separately by each transducer element 12 
and a sample of the magnitude of the echo signal at a 
particular point in time represents the amount of reflec- 
tion occurring at a specific range. Due to differences in 



the propagation paths between a reflecting point and 
each transducer element 12, the echo signals will not 
be detected simultaneously and their amplitudes will not 
be equal. Receiver 16 amplifies the separate echo sig- 

5 nals via a respective TGC amplifier 28 in each receive 
channel. The amplified echo signals are then fed to a 
receive beamformer 30. Each receiver channel of the 
receive beamformer is coupled to a respective one of 
transducer elements 12 by a respective TGC amplifier 

10 28. 

[0027] Under the direction of master controller 20, re- 
ceive beamformer 30 tracks the direction of the trans- 
mitted beam, sampling the echo signals at a succession 
of ranges along each beam. Receive beamformer 30 im- 

15 parts the proper time delay to each amplified echo sig- 
nal, provides dynamic apodization on receive and sums 
the delayed and apodized echo signals to provide a 
summed echo signal which accurately indicates the total 
ultrasonic energy reflected from a point located at a par- 

20 ticular range along one ultrasonic beam. The receive fo- 
cus time delays are computed in real-time using spe- 
cialized hardware or are read from a look-up table. The 
receive channels also have circuitry for filtering the re- 
ceived pulses. The time-delayed receive signals are 

2S then summed and supplied to a signal processor 32. 
Signal processor 32 converts the summed receive sig- 
nals to display data. In the B-mode (grey-scale), this is 
the envelope of the signal with some additional process- 
ing, such as edge enhancement and logarithmic com- 

30 pression. A scan converter 34 receives the display data 
from signal processor 32 and converts the data into the 
desired image for display. In particular, scan converter 
34 converts the acoustic image data from polar coordi- 
nate (R — G) sector format or Cartesian coordinate linear 

35 array to appropriately scaled Cartesian coordinate dis- 
play pixel data at the video rate. These scan-converted 
acoustic data are then provided for display on display 
monitor 22, which images the time-varying amplitude of 
the signal envelope as a grey scale. 

40 [0028] FIG. 2 shows portions of an ultrasound imag- 
ing system in accordance with the present invention. In 
this system each transducer element in the transmit ap- 
erture is pulsed N times by a respective multiphase (e. 
g., bipolar) pulser 24' in accordance with transmit codes 

45 stored in a transmit sequence memory 38. For example, 
the transducer elements are pulsed in accordance with 
a first transmit code during a first transmit firing and in 
accordance with a second transmit code during a sec- 
ond transmit firing, wherein the first and second transmit 

50 codes are applied as phase coding (e.g., polarity revers- 
al) to a conventional transmit pulse. Pulsers 24' drive 
elements 12 of transducer array 10 such that the ultra- 
sonic energy produced is focused at the same transmit 
focal position for each transmit firing. To accomplish this, 

55 identical transmit focus time delays 36 are imparted to 
the respective pulsed waveforms produced by the 
pulsers in accordance with the transmit codes By ap- 
propriately adjusting the transmit focus time delays in a 
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conventional manner, the ultrasonic beams can be fo- 
cused at a multiplicity of transmit foca! positions to effect 
a scan in an image plane. 

[0029] For each transmit, the echo signals from trans- 
ducer elements 12 are fed to respective receive chan- 
nels 40 of the receive beamformer. Under the direction 
of master controller 20 (FIG. 1), the receive beamformer 
tracks the direction of the transmitted beam. The receive 
beamformer imparts the proper receive focus time de- 
lays 42 to the received echo signal and sums them to 
provide an echo signal which accurately indicates the 
total ultrasonic energy reflected from a particular trans- 
mit focal position along a transmit beam. The time-de- 
layed receive signals are summed in a receive summer 
44 for each of the N transmit firings focused at a partic- 
ular transmit focal position. The summed receive signal 
for each of the N transmit firings is then provided in suc- 
cession to a "slow-time" filter 46, which filters across the 
N transmit firings and then supplies a filtered signal to 
signal processor 32. Signal processor 32 forms the en- 
velope of this filtered signal. After post-processing (in- 
cluding edge enhancement and logarithmic compres- 
sion) and scan conversion, a scan line is displayed on 
display monitor 22 (FIG. 1). This procedure is repeated 
so that a respective scan line is displayed for each trans- 
mit focal position (in the case of one transmit focal po- 
sition for each beam angle) or for each transmit vector 
(in the case of multiple transmit focal positions for each 
transmit vector). 

[0030] In accordance with the preferred embodiments 
of the invention, "slow-time" filter 46 comprises an FIR 
filter 48 having an input coupled to the output of receive 
summer 44, and a vector summer 50 having an input 
coupled to FIR filter 48 and an output coupled to signal 
processor 32. The FIR filter has M filter taps for receipt 
of a respective set of M filter coefficients for each trans- 
mit firing. The filter coefficients for the n-th transmit firing 
are a n c Qt a n c^, a n c M _^ , where a n is the scalar weight- 
ing for the n-ih transmit firing, n = 0, 1 , N- 1 , and c 0 , 
c v ■••« c m- 1 ' s a se * °f f'K er coefficients which are select- 
ed so that FIR filter 48 passes a desired frequency band 
in the receive signal. The scalar weightings a 0 , a v .... 
a N .^ t cause the "slow-time" filter to selectively pass or 
attenuate the bandpassed signals as a function of the 
harmonic mode and the scatterer velocities. The filter 
coefficients a n c 0% a n c^, .... a^^ are supplied to the fil- 
ter for each transmit firing by the master controller from 
a filter coefficient memory 52. For example, for the first 
transmit firing, the set of filter coefficients a 0 Co, apC,, .... 
a o c A/M iS supplied to the FIR filter; for the second trans- 
mit firing, the set of filter coefficients a 1 c 0 , a^c v 
a i c M-i is supplied to the FIR filter; and so forth. The filter 
coefficients are programmable depending upon the di- 
agnostic application. Different sets of filter coefficients 
can be stored in look-up tables inside the master con- 
troller memory and the desired set of coefficients can be 
selectable by the system operator. For applications 
where the number of transmit firings N- 2, one or more 



sets of filter coefficients are stored in memory, one set 
of filter coefficients being transferred to the FIR filter be- 
fore the first transmit firing and another set of filter co- 
efficients being transferred to the FIR filter after the first 
5 transmit firing and before the second transmit firing 
(when the same scalar weighting applies to two transmit 
firings, the same filter coefficient set can be used for 
both firings). Similarly, for applications where the 
number of transmit firings N> 2, two or more sets of filter 
10 coefficients are stored in memory. The successive FIR 
filter output signals for the N transmit firings are accu- 
mulated in a vector summer 50. The output signal of the 
vector summer then undergoes conventional B-mode 
processing, followed by scan conversion and display. 
15 [0031] In accordance with a preferred embodiment of 
the invention, firing-to-firing (i.e., "slow-time") filtering is 
combined with transmit phase coding to produce an en- 
hanced ultrasound image. The "slow-time" filter re- 
sponds differently to the different modes (fundamental, 
20 second subharmonic, second harmonic, third harmonic, 
etc. ) of the reflected signal because the transmit phases 
change over the set of transmit firings. This permits one 
to design the transmit phases and "slow-time" filter to 
selectively enhance the desired modes while suppress- 
es ing others. In particular, if the transmitted signal has a 
phase term expt/GJ, where /= 0, 1 , .... A/- 1 , then the k- 
th (sub)harmonic has a phase term exp[y/c<- 1 >G J, / - 0, 
1 , . . ., N - 1 . Thus, if the "slow-time" filter coefficients are 
/ = 0, 1 , .... N- 1, then the effective "slow-time" filter 
30 for the /c-th (sub)harmonic is apxp\jk(^ )ft J, / = 0, 1 , 
N - 1 , which has a transfer function that depends on the 
particular mode k. 

[0032] The system of the present invention has three 
different application areas; contrast harmonic imaging, 

35 tissue harmonic imaging and B-flow imaging. For each 
application, the transmit phases and "slow-time" filter 
weightings may be selected to achieve the desired fil- 
tering on the fundamental and the (sub)harmonics. The 
"slow-time" filter responses for various preferred em- 

40 bodiments are shown in FIGS. 4-1 1 . The "slow-time" fil- 
ter response to the fundamental mode is indicated by 
solid lines, to the second harmonic is indicated by 
dashed lines, and to the second subharmonic is indicat- 
ed by dotted lines. The horizontal axis corresponds to 

45 "slow-time" normalized frequency, while the vertical axis 
is the magnitude of the "slow-time" filter output. The ex- 
pected range of normalized operating frequencies lies 
in the range of ±0.2. 

[0033] In contrast harmonic imaging, contrast agents 
50 made up of gas-filled microbubbles are injected into the 
blood to serve as markers for imaging blood flow. It is 
desired to see fundamental or second harmonic signals 
from contrast flow with suppressed background tissue 
signals and little motion flash artifacts. The broad-band 
55 pulses transmitted to a particular transmit focal position 
in sequence are phase coded. In particular, N pulses 
centered at frequency f 0 are transmitted to each trans- 
mit focal position. On receive, a "slow-time" filter ex- 
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tracts the (sub)harmonic flow signal over the N trans- 
mits. In particular, a set of n slow-time° filter weightings 
a 0 , a,, .... a M _ v are selected so that the /W-tap "slow- 
time 0 FIR filter 48 passes substantially all of the desired 
harmonic or subharmonic frequencies in signals reflect- s 
ed from agents moving at certain velocities, while sub- 
stantially suppressing signals at the fundamental fre- 
quencies. If the transmitted center frequency is at 
then tissue/contrast nonlinearities will generate har- 
monics at k/ Q , where k is an integer greater than or io 
equal to 2. Also, subharmonics at frequencies /<//cmay 
be generated by contrast bubble destruction. 
[0034] In accordance with one preferred embodiment 
of the invention, contrast harmonic imaging is achieved 
by high-pass filtering the second harmonic and sup- is 
pressing all of the fundamental signal, which results in 
good suppression of background tissue signals, as seen 
in FIG. 6. The response shown in FIG. 6 was obtained 
using transmit phases [0°, 180°, 180°, 0°] and filter 
weightings [0.4, 1,-1, -0.4]. FIG. 6 shows the case 20 
wherein, in addition to suppression of the fundamental, 
stationary components of the (sub)harmonic signal are 
substantially blocked by a notch filter. This preferred em- 
bodiment is useful for nonlinear imaging of nonstation- 
ary or flow regions, and in particular for imaging contrast 2s 
agents injected into the blood stream. 
[0035] In another preferred embodiment of the inven- 
tion, contrast harmonic imaging is realized by high-pass 
filtering of both the fundamental and the second har- 
monic signals, which results in better low-flow sensitivity 30 
but larger motion flash artifacts, as seen in FIG. 11 . The 
response shown in FIG. 11 was obtained using transmit 
phases [0°, 180°, 0°, 180°] and filter weightings [1.1, 
-1.-1]. 

[0036] In other preferred embodiments of the inven- 3S 
tion, contrast harmonic imaging is realized by high-pass 
filtering or suppressing the fundamental signal and all- 
pass filtering the second harmonic signal, which results 
in more tissue background (from the second harmonic) 
but shows harmonic signals from even the slowest-mov- 40 
ing contrast agents. Examples of this mode are shown 
in FIGS. 4, 5 and 7-10. The response shown in FIG. 4 
was obtained using transmit phases [0°, 180°, 0°, 180°] 
and filter weightings [0.4, 1,1, 0.4]; the response shown 
in FIG. 5 was obtained using transmit phases [0°, 90°, 45 
0°, 180°], filter weightings (0.4, 1,1, 0.4] and filter phas- 
es [0°, 90° : 0°, 0°]; the response shown in FIG. 7 was 
obtained using transmit phases [0°, 180°] and filter 
weightings [1 , 1 ]; the response shown in FIG. 8 was ob- 
tained using transmit phases [180°, 0°, 180°] and filter so 
weightings [0.5, 1, 0.5]: the response shown in FIG. 9 
was obtained using transmit phases [0°, 0°, 180°, 180°] 
and filter weightings [1 , 1 , 1 , 1 ]; and the response shown 
in FIG. 10 was obtained using transmit phases [0°, 180°, 
180° : 0°] and filter weightings [1, 1,1, 1]. The response S5 
shown in FIG. 5 was obtained using a complex filter. 
[0037] In tissue harmonic imaging, the goal is to see 
harmonic signals (in particular, the second harmonic) 



generated by nonlinear propagation in tissue. In yet an- 
other preferred embodiment of the invention, this is 
achieved by suppressing all of the fundamental signal 
and passing all of the second harmonic. To this end, the 
transmit phase codes and "slow-time" filter weightings 
represented by FIGS. 4, 5, 8 or 10 can be used. The 
transmit phase codes and "slow-time" filter weightings 
which produce the response shown in FIG. 7 can also 
be used, but with larger flash motion artifacts. 
[0038] Finally, the goal in B-mode flow imaging is to 
visualize fundamental signals from blood flow (without 
contrast agents) with minimal motion flash artifacts. In 
a preferred embodiment of the invention, this is accom- 
plished by high-pass filtering the fundamental and all- 
pass filtering the second harmonic. Passing the second 
harmonic smoothes out the flash artifacts. To this end, 
the transmit phase codes and "slow-time" filter weight- 
ings represented by FIGS. 7 or 9 can be used. 
[0039] In accordance with a further preferred embod- 
iment of the invention, the B-mode flow image is super- 
imposed on a conventional B-mode image. This allows 
the diagnostician to observe the flow of blood relative to 
known anatomical landmarks during medical diagnosis. 
This B-mode image feed-through is achieved by per- 
turbing one of the "slow-time" filter weightings. For ex- 
ample, the weighting Sq for the first transmit firing (or for 
any other transmit firing) can be perturbed by an amount 
a, as shown in FIG. 3. The B-mode feed-through allows 
the flow image to be superimposed on top of a conven- 
tional B-mode image for display. Alternatively, the flow 
image may be superimposed in color on top of a con- 
ventional B-mode image for display. 
[0040] The time interval between each of the N trans- 
mits per focal position is user controllable to determine 
the "slow-time" filter cut-off frequency. A longer interval 
between each of the A/transmits to a particular focal po- 
sition results in a lower cutoff frequency with higher sen- 
sitivity to low velocity flow. 



Claims 

1. A system for imaging ultrasound scatterers, com- 
prising: 

an ultrasound transducer array comprising a 
plurality of transducer elements; 
pulsing means coupled to said transducer array 
for pulsing selected ones of said transducer el- 
ements which form a transmit aperture with a 
first phase-coded transmit pulse as a function 
of a first transmit phase code during said first 
transmit firing and with a second phase-coded 
transmit pulse as a function of a second trans- 
mit phase code during said second transmit fir- 
ing; 

transmit beamforming means coupled to said 
pulsing means for forming first and second 
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beams during said first and second transmit fir- 
ings respectively, said first and second beams 
being focused at substantially the same trans- 
mit focal position; 

receive beamforming means coupled to said 5 
transducer array for forming a first beam- 
summed receive signal from a first set of re- 
ceive signals from other selected ones of said 
transducer elements which form a receive ap- 
erture subsequent to said first transmit firing io 
and forming a second beamsummed receive 
signal from a second set of receive signals from 
said other selected ones of said transducer el- 
ements which form said receive aperture sub- 
sequent to said second transmit firing; 15 
filtering means for forming a first filtered signal 
by applying a first "slow-time" filter weighting a Q 
to said first beamsummed receive signal and 
for forming a second filtered signal by applying 
a second "slow-time' 1 filter weighting at to said 20 
second beamsummed receive signal; 
a vector summer coupled to said filtering 
means for summing at least said first and sec- 
ond filtered signals to form a "slow-time" filtered 
signal; 2s 
means for processing said "slow-time" filtered 
signal to form an image signal; and 
means for displaying an image which is a func- 
tion of said image signal. 

30 

2. The system of claim 1 wherein said filtering means 
comprise: 

means for supplying first and second sets of fil- 
ter coefficients, said first set of filter coefficients 35 
being derived by multiplying a predetermined 
set of coefficients by said first "slow-time" filter 
weighting a 0 and said second set of filter coef- 
ficients being derived by multiplying said pre- 
determined set of coefficients by said second 40 
"slow-time" filter weighting ay, and 
a filter having a signal input coupled to an out- 
put of said receive beamforming means, a plu- 
rality of filter taps coupled to receive said first 
and second sets of filter coefficients from said 45 
means for supplying first and second sets of fil- 
ter coefficients, and an output for supplying said 
first filtered signal in dependence on said first 
beamsummed receive signal and said first set 
of filter coefficients and for supplying said sec- so 
ond filtered signal in dependence on said sec- 
ond beamsummed receive signal and said sec- 
ond set of filter coefficients. 

3. The system of claim 1 wherein a 0 = -a x . ss 

4. The system of claim 1 wherein said pulsing means 
are adapted to pulse said selected transducer ele- 



ments which form said transmit aperture with a third 
phase-coded transmit pulse as a function of a third 
transmit phase code during a third transmit firing; 
said transmit beamforming means being adapted to 
form a third beam during said third transmit firing, 
said third beam being focused at said transmit focal 
position; said receive beamforming means being 
adapted to form a third beamsummed receive signal 
from a third set of receive signals from said selected 
transducer elements which form said receive aper- 
ture subsequent to said third transmit firing; said fil- 
tering means being adapted to form a third filtered 
signal by applying a third "slow-time" filter weighting 
to said third beamsummed receive signal; and said 
vector summer being adapted to sum at least said 
first, second and third filtered signals to form said 
"slow-time" filtered signal. 

5. The system of claim 4 wherein said first through 
third transmit phase codes are [180°, 0°, 180°] and 
said first through third "slow-time" filter weightings 
are [0.5, 1, 0.5]. 

6. The system of claim 4 wherein said pulsing means 
are adapted to pulse said selected transducer ele- 
ments which form said transmit aperture with a 
fourth phase-coded transmit pulse as a function of 
a fourth transmit phase code during a fourth trans- 
mit firing; said transmit beamforming means are 
adapted to form a fourth beam during said fourth 
transmit firing, said fourth beam being focused at 
said transmit focal position; said receive beamform- 
ing means are adapted to form a fourth beam- 
summed receive signal from a fourth set of receive 
signals from said selected transducer elements 
which form said receive aperture subsequent to 
said fourth transmit firing; said filtering means are 
adapted to form a fourth filtered signal by applying 
a fourth "slow-time" filter weighting to said fourth 
beamsummed receive signal; and said vector sum- 
mer is adapted to sum at least said first through 
fourth filtered signals to form said "slow-time" fil- 
tered signal. 

7. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 180°, 0°, 180°] 
and said first through fourth "slow-time" filter 
weightings are [0.4, 1,1, 0.4]. 

8. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 90°, 0°, 180°], 
said first through fourth "slow-time" filter weightings 
are [0.4, 1,1, 0.4], and first through fourth filter 
phases are [0°, 90°, 0°, 0°], said first through fourth 
filter phases being respectively adapted to be ap- 
plied to said filtering means in conjunction with said 
first through fourth "slow-time" filter weightings. 
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9. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 180°, 180°, 0°] 
and said first through fourth "slow-time" filter 
weightings are [0.4, 1,-1, -0.4]. 

10. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 0°, 180°, 180°] 
and said first through fourth "slow-time" filter 
weightings are[1, 1, 1, 1]. 

11. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 180°, 180°, 0°] 
and said first through fourth "slow-time" filter 
weightings are[1, 1, 1, 1]. 

12. The system of claim 6 wherein said first through 
fourth transmit phase codes are [0°, 1 80°, 0°, 1 80°] 
and said first through fourth "slow-time" filter 
weightings are [1 , 1 , -1 , -1 ]. 

13. The system of claim 6 wherein said first through 
fourth transmit phase codes and said first through 
fourth "slow-time" filter weightings are selected so 
that said filtering means can high-pass filter a sec- 
ond harmonic signal and substantially suppress a 
fundamental signal. 

14. The system of claim 6 wherein said first through 
fourth transmit phase codes and said first through 
fourth "slow-time" filter weightings are selected so 
that said filtering means can high-pass filter a fun- 
damental signal and a second harmonic signal. 

15. The system of claim 6 wherein said first through 
fourth transmit phase codes and said first through 
fourth "slow-time" filter weightings are selected so 
that said filtering means can all-pass filter a second 
harmonic signal and substantially suppress a fun- 
damental signal. 

16. The system of claim 6 wherein said first through 
fourth transmit phase codes and said first through 
fourth "slow-time" filter weightings are selected so 
that said filtering means can all-pass filter a second 
harmonic signal and high-pass filter a fundamental 
signal. 

17. A method for imaging ultrasound scatterers, com- 
prising the steps of: 

selecting a first transmit phase code and a first 
"slow-time" filter weighting for a first transmit fir- 
ing, and a second transmit phase code and a 
second "slow-time" filt r weighting for a second 
transmit firing; 

driving a first set of transducer elements form- 
ing a transmit aperture in a transducer array 
with a first phase-coded transmit pulse as a 



function of said first transmit phase code during 
said first transmit firing, said first transmit beam 
being focused at a transmit focal position; 
receiving a first set of echo signals from a sec- 

s ond set of transducer elements forming a re- 

ceive aperture in the transducer array subse- 
quent to said first transmit firing; 
forming a first beamsummed receive signal 
from said first set of echo signals; 

io forming a first filtered signal by applying a first 

"slow-time" filter weighting to said first beam- 
summed receive signal; 
driving said first set of transducer elements with 
a second phase-coded transmit pulse as a 

15 function of said second transmit phase code 

during said second transmit firing, said second 
transmit beam being focused at said transmit 
focal position; 

receiving a second set of echo signals from said 
20 second set of transducer elements subsequent 

to said second transmit firing; 

forming a second beamsummed receive signal 

from said second set of echo signals; 

forming a second filtered signal by applying a 
25 second "slow-time" filter weighting a 1 to said 

second beamsummed receive signal; 

summing at least said first and second filtered 

signals to form a "slow-time" filtered signal; 

processing said "slow-time" filtered signal to 
30 form an image signal; and 

displaying an image which is a function of said 

image signal. 

18. The method of claim 17 wherein said first filtered 
35 signal is a function of said first beamsummed re- 
ceive signal and a first set of filter coefficients de- 
rived by multiplying a predetermined set of coeffi- 
cients by said first "slow-time" filter weighting a^, 
and said second filtered signal is a function of said 
40 second beamsummed receive signal and a second 
set of filter coefficients derived by multiplying said 
predetermined set of coefficients by said second 
"slow-time" filter weighting a v 

45 19. The method of claim 17 wherein a 0 = -a v 

20. A method for imaging ultrasound scatterers in a 
flowing fluid medium, comprising the steps of: 

50 selecting a first transmit phase code and a first 

"slow-time" filter weighting for a first transmit fir- 
ing, and a second transmit phase code and a 
second "slow-time" filter weighting for a second 
transmit firing; 

55 injecting a contrast agent into the flowing fluid 

medium; 

driving a first set of transducer elements form- 
ing a transmit aperture in a transducer array 
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with a first phase-coded transmit pulse as a 
(unction of said first transmit phase code during 
satd first transmit firing, said first transmit beam 
being focused at a transmit focal position: 
receiving a first set ol echo signals from a sec- 5 
ond set of transducer elements forming a re- 
ceive aperture in the transducer array subse- 
quent to said first transmit firing; 
forming a lirst beamsummed receive signal 
from said first set of echo signals; 10 
forming a first filtered signal by applying a first 
•slow-time" filter weighting a 0 to said first beam- 
summed receive signal; 

driving said first set of transducer elements with 
a second phase-coded transmit pulse as a is 
function of said second transmit phase code 
during said second transmit firing, said second 
transmit beam being focused at said transmit 
focal position 

receiving a second set of echo signals from said 20 

second sot ol ttansducei elements subsequent 

to said second transmit firing 

forming a second ocamsummed receive signal 

from said second set of echo signals; 

forming n seccod ftrtcrcd signal by applying a 2s 

second "slcw-timo* filtor weighting a 1 to said 

second bcamsjmmed receive signal; 

summing at tonst said trst and second filtered 

signals to form a 'slow-time* filtered signal; 

processing said •slow-time" filtered signal to 30 

form an imago sipjnal and 

displaying an imago when is a function of said 

image signal 

21. A system for imaging ultrasound scatterers, com- 35 
prising: 

a transducer a'ray comprising a plurality of ul- 
trasound transducer elements: 
means for conirollmg said transducer array to *o 
transmit a first transmit ultrasound beam having 
a first phase toward said ultrasound scatterers 
during a first transmit firing and to transmit a 
second transmit ultrasound beam having a sec- 
ond phase different than said first phase toward 45 
said ultrasound scatterers during a second 
transmit firing: 

means for acquiring data from said transducer 
array representing a fiist receive ultrasound 
beam derived from reflections of said first trans- so 
mit ultrasound beam from said ultrasound scat- 
terers and data from said transducer array rep- 
resenting a second receive ultra-sound beam 
derived from reflections of said second transmit 
ultrasound beam from said ultrasound scatter- 55 
ers; 

means for filtering said first receive ultrasound 
beam data to form a first filtered signal following 



said first transmit firing and for filtering said sec- 
ond receive ultrasound beam data to form a 
second filtered signal following said second 
transmit firing, said first filtered signal being a 
function of a first "slow-time" filter weighting 
and said second filtered signal being a function 
of a second "slow-time" filter weighting; 
means for summing at least said first and sec- 
ond filtered signals to form a "slow-time" filtered 
signal; 

processing means for producing an image sig- 
nal which is a function of said "slow-time" fil- 
tered signal; and 

means for displaying an image which is a func- 
tion of said image signal. 

22. The system of claim 21 wherein: 

said transducer array controlling means is 
adapted to control said transducer array to 
transmit a third transmit ultrasound beam hav- 
ing a third phase toward said ultrasound scat- 
terers during a third transmit firing; 
said data acquiring means is adapted to ac- 
quire data from said transducer array repre- 
senting a third receive ultrasound beam derived 
from reflections of said third transmit ultra- 
sound beam from said ultrasound scatterers; 
said filtering means is adapted to filter said third 
receive ultrasound beam data to form a third 
filtered signal following said third transmit firing, 
said third filtered signal being a function of a 
third "slow-time" filter weighting; and 
said summing means is adapted to sum at least 
said first, second and third filtered signals to 
form said "slow-time" filtered signal. 

23. The system of claim 22 wherein: 

said transducer array controlling means is 
adapted to control said transducer array to 
transmit a fourth transmit ultrasound beam hav- 
ing a fourth phase toward said ultrasound scat- 
terers during a fourth transmit tiring; 
said data acquiring means is adapted to ac- 
quire data from said transducer array repre- 
senting a fourth receive ultra-sound beam de- 
rived from reflections of said fourth transmit ul- 
trasound beam from said ultrasound scatterers; 
said filtering means is adapted to filter said 
fourth receive ultrasound beam data to form a 
fourth filtered signal following said fourth trans- 
mit firing, said fourth filtered signal being a func- 
tion of a fourth "slow-time" filter weighting: and 
said summing means is adapted to sum at least 
said first through fourth filtered signals to form 
said "slow-time" filtered signal. 
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